Classification of general n-qubit states under stochastic local operations and classical 
communication in terms of the rank of coefficient matrix 
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We solve the entanglement classification under stochastic local operations and classical commu- 
nication (SLOCC) for general n-qubit states. For two arbitrary pure n-qubit states connected via 
local operations, we establish an equation between the two coefficient matrices associated with the 
states. The rank of the coefficient matrix is preserved under SLOCC and gives rise to a simple 
way of partitioning all the pure states of n qubits into different families of entanglement classes, as 
exemplified here. When applied to the symmetric states, this approach reveals that all the Dicke 
states \£, n) with £ = !,..., [n/2] are inequivalent under SLOCC. 



I. INTRODUCTION 

Entanglement, a key feature that distinguishes quan- 
tum information from classical information, has applica- 
tions in cryptography, teleportation and quantum com- 
putation [lj. While the bipartite entanglement is well 
understood, the task of classifying multipartite entangle- 
ment beyond two qubits becomes increasingly difficult. 
To classify entangled states, some equivalence relation 
has to be introduced. Of particular importance is the 
equivalence under stochastic local operations and classi- 
cal communication (SLOCC) since two states belonging 
to the same equivalence class can perform the same tasks 
of quantum information theory. 

For three qubits, in terms of the local ranks of the 
reduced density matrices, it has been shown that there 
are six inequivalent SLOCC classes Q. For four or more 
qubits, there exists an infinite number of inequivalent 
SLOCC classes. It is highly desirable to partition the in- 
finite SLOCC classes into a finite number of families such 
that states belonging to the same family possess simi- 
lar properties, according to some criteria for determining 
which family a given state belongs to. Considerable ef- 
forts have been undertaken over the last decade for the 
SLOCC entanglement classification of four-qubit states 
resulting in a finite number of families [3-6] or classes 
0-[Hj]- For more than four qubits, a few attempts have 
been made for SLOCC classification for subsets of the 
general n-qubit states such as the Greenberger-Horne- 
ZeilingerJGHZ)-type, W-type, and GHZ-W-type n-qubit 



states [13j, symmetric n-qubit states [14H16B . even n- 
qubit states [T3, , and odd n-qubit states (13] ■ Despite 
these efforts, a SLOCC classification for general n-qubit 
states is still beyond reach. 

Our aim is to solve the problem of SLOCC classifica- 
tion for all multipartite pure states in the general n-qubit 
case. To this end, we demonstrate that the rank of the co- 
efficient matrix of a pure n-qubit state is invariant under 
SLOCC. SLOCC invariants for subsets of n-qubit states 



have been the subject of several recent studies [17h21|. 
In 17|, LL8| , the invariant element is the determinant of 
coefficient matrices of even n qubits. In (l9^ . the invari- 
ant element is the rank of square matrices of size two 
constructed using three functions defined on the space of 
odd n qubits. 

We construct the coefficient matrices for general n- 
qubit states by arranging the coefficients in lexicograph- 
ical order. For two states connected via local operations, 
their coefficient matrices are related through an equa- 
tion. In the case where the local operations are invert- 
ible, the two states are said to be SLOCC equivalent and 
the two coefficient matrices have the same rank; i.e., the 
rank is preserved under SLOCC. The rank gives rise to 
a simple way of partitioning all the n-qubit states into 
different SLOCC families. For n-qubit symmetric Dicke 
states \£, n) with I (I = 1, . . . , [n/2]) excitations, we show 
that the rank of the coefficient matrix of \£, n) is equal to 
£+1 and, therefore, all these states are inequivalent under 
SLOCC. Finally, composing the rank and permutations 
of qubits allows us to define subfamilies by cutting each 
family in pieces. 



II. SLOCC MATRIX EQUATION AND THE 
INVARIANCE OF THE RANK 



We write the state of n qubits as \tp') = 

^^ =0 1 ai|i), where |z) are basis states and a>i are co- 
efficients. We associate to an n-qubit state a 
2 [n/2] x 2 [(«+i)/2] coefficient matrix M{\yb')) whose en- 
tries are the coefficients ao, a±, . . . , a 2 ^~\ arranged in as- 
cending lexicographical order. To illustrate, we list below 
M(\ip')) for n = 3: 



M{W)) 



a ai a 2 a 3 
04 as a,Q a-j 



(1) 
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and for n = 4: 



M(|V0) = 



^ a oi a 2 a 3 ^ 

a4 0,5 a6 a7 

as ag aio an 

\ 012 Oi 3 O14 015 / 



(2) 



We refer to the rank of the coefficient matrix M(\ip'}) 
as the rank of the state denoted as rank(|V> ))• 

We exemplify with the n-qubit |GHZ) state ^(|0)® n + 
and we find rank(|GHZ)) = 2. It is clear that the 
rank of any n-qubit state ranges over the values 1, 2, . . . , 

2 [n/2] < 

Theorem. Let \ip) be another state of n qubits with 

\ip) = J2i=o 1 an( i ^(IV*)) be the corresponding co- 
efficient matrix constructed in the same manner as was 
done for M(\tp')). If the states \tp) and \tp') are related 
by 



\ip') =Ai®A 2 ®-.- ®A n \i)), 



(3) 



where the local operators A\ , A2 , ■ • • , and „4„ are not 
necessarily invertible, then the following matrix equation 
holds for general n qubits: 

M(|V'» = (At®. ■ ■ ®-4[»/2])M(|V))(-4[„/2]+i®- ■ ■ ®-4n) T 

(4) 

Equation JU) holds particularly true for two SLOCC 
equivalent states \t/j) and \4>') which satisfy Eq. ^ along 
with the local operators A\, A2 , ■ • ■ , and A n being in- 
vertible [2]. It follows from Eq. gj) that two SLOCC 
equivalent states have the same rank, in other words, the 
rank is invariant under SLOCC, thereby revealing that 
the rank is an inherent physical property. Therefore, if 
two states differ in their ranks, then they belong neces- 
sarily to different SLOCC equivalent classes. It should 
be noted that the converse does not hold; i.e., two states 
with the same rank are not necessarily equivalent. 

We have the following two simple results: (i) The rank 
of a full separable state is always 1. (ii) The rank of a 
genuinely entangled state is always greater than 1. 

Remark 1. Taking the determinants of both sides of 
Eq. ((4]) for even n yields 

[13, El: 



detM(|V>')) = det M(|V))[det(^i) . . . det(X)] 2 

(5) 

It follows from Eq. © that if one of detMfl^')) and 
det M( I ■)/>)) vanishes while the other does not, then the 
state is not equivalent to under SLOCC. In view 
of the fact that the determinant of a matrix is nonvanish- 
ing if and only if it has full rank, the SLOCC invariance 
of the rank is stronger than the invariance of the deter- 
minant. 



III. SLOCC CLASSIFICATION IN TERMS OF 
THE RANK 

We define the family T n>r to be the set of all n-qubit 
states with the same rank r. In the sequel, we will omit 
the subscript n and simply write J>, whenever the num- 
ber of qubits is clear from the context. Thus, there exist 
2[™/ 2 I different SLOCC families for any n qubits. Clearly, 
if two states are SLOCC equivalent then they belong to 
the same family. However, the converse does not hold: 
two states belonging to the same family may be inequiv- 
alent under SLOCC. It is further noted that when n > 4, 
at least one family contains an infinite number of SLOCC 
classes. 

For any n qubits, the following hold: (i) Family T\ 
contains all the full separable states, (ii) Family T\ con- 
tains no genuine entangled states, (iii) Family T\ con- 
tains finite SLOCC classes, (iv) Family F2 contains the 
n-qubit |GHZ) state, (v) Family J 7 2+r contains the fol- 
lowing state: 

-J=(|0) - |2« - 1} + £ |fc(2 [(n+1)/2] + 1))), (6) 

where 1 < r < 2^/^ - 2. 

Now we turn to the n-q ubit symmetric Dicke states 
\£, n) with £ excitations [221 ]: 

% n) = QT 1/2 E Pk\U, la, ... , U, O1+1, ■ • ■ , 0»>, (7) 



where I ranges from 1 to n — 1 and {Pk} is the set of 
all distinct permutations of the spins. These states have 
been featured in theoretical studies 0, [24| and imple- 
mented experimentally [25|, [26|. The Dicke state |l,n) 
is just the n-qubit \W) state and \£,n) is equivalent to 
\n-i,n) under SLOCC. 

As has been previously noted, all symmetric Dicke 
states \£,n) with I = l,...,[n/2] are SLOCC inequiv- 



i) wi 

m 



alent [LJ, Il5j . These states, as demonstrated below, 
can also be distinguished by the rank of their coefficient 
matrices which depends only on the number of excita- 
tions and is independent of the number of qubits. Since 
rank(|£, n)) =rank(|n — £, n)), we only need to compute 
rank(|^,n)) with 1 < I < [n/2]. This can be done as fol- 
lows. We first construct the coefficient matrix M(\ip')) 
of state in the same manner as discussed above. 
We may write \if>) in terms of an orthogonal basis as 
|V/) = J2 a tn2---i n 1*1*2 •■ -in), where iii 2 ...i„ is the 71- 
bit binary form of the index i. Inspection of the struc- 
ture of the matrix M(\ip')) reveals that the coefficient 



a 



n— *r„/ai*[n/ai+i— *' 



is the entry in the (ii ■ ■ ■ ir„/2])th row 
and (£[71/21+1 • • • £n)th column of the matrix. Here, the n 
bits are split into two halves, referred to as the row bits 
and column bits, respectively: bits 1 to [n/2] are used 
to specify the row number, and bits [n/2] + 1 to n are 
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used to specify the column number. In view of Eq. (J7J, 
the nonzero entries of the coefficient matrix M(\£, n)) are 
those whose n-bit binary forms have £ bits equal to 1 and 
the rest of the bits equal to 0. We further observe that the 
rows of M(\£, n)) with no more than £ row bits equal to 
1 are nonzero rows, while the remaining rows are identi- 
cally zero. Consider the rows with j (0 < j < £) row bits 
equal to 1. Clearly, there are ("^ 2 ) such rows of M(\£, n)) 
that are identical. Letting j vary from to £ gives a total 
of £ + 1 different rows. It can be verified that these £+1 
rows are independent. This yields rank(\£, n)) = £ + 1. 

Accordingly, for any n qubits, all the Dicke states \£, n) 
with £ = 1, [n/2], are inequivalent under SLOCC, 
since they differ in their ranks. Further, it can readily 
be seen that the Dicke state \£, n) with £ = 1, . . . , [n/2] 
belongs to the family J-£+i- This gives rise to a complete 
SLOCC classification of all the symmetric Dicke states. 

Remark 2. It follows from the discussion above that 
both the states \W) and |GHZ) have the same rank, 
thereby revealing that the two states have a similar alge- 
braic structure. It is further noted that both the states 
\W) and | GHZ) admit a similar Frobenius algebra struc- 
ture m. 



IV. RANKS OF COEFFICIENT MATRICES 
UNDER PERMUTATIONS OF QUBITS 



In [Tfl, [18] , we presented a systematic way to find all 
the possible coefficient matrices for even n-qubit states 
such that the determinants of these coefficient matrices 
are invariant under SLOCC. Here we extend this con- 
struction to general n qubits. Observe that to write a 
2[n/2] x 2lO+i)/2] ma trix into binary index form, we need 
[n/2] row bits and [(n+ l)/2] column bits. In the binary 
form of the coefficient matrix given in Eqs. (Q} and ([2]), 
bits 1 to [n/2] are the row bits, and bits [n/2] + 1 to n are 
the column bits. Alternatively, we may choose any [n/2] 
bits as the row bits and the remaining [(n + l)/2] bits 
as the column bits. This amounts to (±) n+1 mod 2 ( [n ™ 2] ) 
different coefficient matrices, ignoring those that end up 
exchanging rows or columns. The factor of 1/2 for even 
n arises because exchanging the row and column bits is 
equivalent to transposing the matrix. It turns out that 
the ranks of these coefficient matrices are all invariant 
under SLOCC. To see this, we will resort to permuta- 
tions of qubits. Let a be a permutation of qubits given 
by Q 



o- = (qi,t 1 )(q 2 ,t 2 ) ■ ■ ■ (qk,t k ), 



(8) 



where {qi,U) is the transposition of a pair of qubits qi 
and U with qt being a row bit and t% a column bit. Ex- 
hausting all possible values of qi, . . . , qu and t\ , . . . , tk 
such that 1 < q\ < q 2 < ... < q k < [(n + l)/2], 
[n/2] < ti < <2 < . . . < tk < n, and letting k vary 



from to [{n — l)/2] (we define a = I for k = 0), yields 
7) 2 ,) different permutations of qubits. Let 



^ l^n+l mod 2 | 



M a (\ip'}) denote the coefficient matrix of the state 
under permutation cr, and let ranked?//)) denote its rank. 
We may refer to ranked?/ 1 ')) as the rank of the state \ip'} 
under permutation a. Simply taking the permutation 
a to both sides of Eq. (gj) yields the following SLOCC 
matrix equation: 

M ff (|V/)) - (Ax(i) ® • ■ • ® A a{[n/2]) )M a {m 

(Ar([n/2] + l) ® • • • ® At(«)) T - (9) 

It follows immediately from Eq. ([9]) that two SLOCC 
equivalent states have the same rank with respect to the 
permutation a. That is, the rank with respect to the per- 
mutation a is invariant under SLOCC. Conversely, if two 
states differ in their ranks with respect to the permuta- 
tion cr, then they belong necessarily to different SLOCC 
classes. We define the family T% to be the set of all n- 
qubit states with the same rank r with respect to the 
permutation cr, where r ranges from 1 to 2 <vn / 2 ^ and we 
have omitted a subscript n. Suppose ai,a 2 , ■ ■ ■ , &m with 
m < {\) n+1 mod 2 ([n/2]) ' s a secmence °f permutations 
of the form given in Eq. (JS]). In terms of the rank 
of M CTl , the n-qubit states are divided into 2^- n / 2 ^ fam- 
ilies: J 7 ^ 1 . Then, in terms of the rank of M a2 , each 
family J 7 ^ 1 can be further divided into 2[™/ 2 l subfami- 
lies: T°l°l = T°l n J7 2 2 . Here, each subfamily is 
the intersection of the families and . Assume 
that in terms of the ranks of M CTl , M° 2 , . . . , M <7m ~ 1 , 
the n-qubit states are divided into 2( m_1 ^"/ 2 l families: 
tfi°Z: a ,r~\ = n ... n KZ-l- Then, in terms 
of the rank of M CT ™, each family Jv^^V.'.Ttw^i can t> e 
further divided into 2l"/ 2 l subfamilies: J 7 ^- a j," m = 
Kl a rl\:Zr~\ n = T°l n . .. n T°™. This gives a 
total of 2 m ["/ 2 ] different SLOCC families. 

We exemplify with the family L a .^, 2 = a(|0000) + 
|1111)) + 6(|0101> + |1010)) + |0011) + |0110) for four 
qubits presented by Verstraete et al. As shown in 
Table [H the family L a2 i> 2 is further divided into four sub- 
families (all other subfamilies are empty) with respect to 
permutations a\ = I and a 2 — (1,4): F^ 1 " 2 , J 7 ^ 1 ^ 2 , and 
J-*42 CT2 contain only a single SLOCC class, while Jjj" 2 
contains an infinite number of SLOCC classes. In a simi- 
lar fashion, we can further divide other families presented 
by Verstraete et al. [3| into subfamilies. 

TABLE I. SLOCC classification of L a2 t 2 



T" 1 



T 17 
J o 







7T D "2 
•^2 



r 3 

ab — & a ^ b 



ab 

T , °'2 
-^4 



TT cr 2 

a = i = 0a = ±()&a/0 a ±b 

-pa ! cr 2 -r<y 1 cr 2 TT a 1 cr 2 T " l ff 2 

•> 2,1 ^3,3 •> 4,2 ^4,3 

a = b = ab = & a^b a = ±b&ia^0ab^0&ia^±b 
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Consider the family span{O fc *, fc *} = |0000) + 
|1100) + a|0011) + /3|1111) for four qubits presented by 
Lamata et al. [5]. As shown in Table UH the family 
{0/ c , L, O/c^} is further divided into four subfamilies (all 
other subfamilies are empty) with respect to permuta- 
tions o-i = I and cr 2 = (1,4): J 7 ?,? 3 , and F^p 
contain only a single SLOCC class, while J 7 ^ 1 ^ 2 contains 
an infinite number of SLOCC classes. In a similar way, 
other families presented by Lamata et al. [f| can also be 
further divided into subfamilies. 



TABLE II. SLOCC classification of span{O fc *, O fc *} 



TTf l 








a = /3 


a^P 








7T (J 2 


7T° - 2 
^2 




7T° - 2 





Q = = 


afi = & a / /3 


a/3 / 


T-CT i GT2 

1,2 


7r CT i (7 2 
'1,4 


1 CT 2 
2,3 


J 2,4 



a = /3 = a = P a,0 = & a / /3 a/3 ^ & a / /3 



By using the filters, it has been shown that four five- 
qubit states |\&2), 1^4), l^s), and are in different 
orbits 28]. Letting oi = /, o 2 = (1,5), and 0-3 = (1,3), 
it can be shown that the above four states belong to 
the families F^f\ re- 

spectively. Therefore, these five-qubit states can also be 
distinguished by ranks. Furthermore, it has been shown 
that five six-qubit states |S 2 ), IS4), |5 5 ), |Sg), and |EV) 
are distinguished by the six-qubit filters [28j . Letting 
<7i =/,o-2 = (1,4), and 0-3 = (1,5), it can be shown 
that the above five states belong to the families J r 2 T 2 ° 2 2 ' T3 , 

2,2,4 > 2,4,4 1 •> 3,4,4 i ancL j- 3 3 3 , respectively. 
Therefore, these six-qubit states can also be distinguished 
by ranks. 



with representative | ( 1 0000) + 1 1 1 1 1 ) - 1 001 1 ) - 1 1 100) ) + 
1 01 10). In light of Theorem 1 in [2(j, the two representa- 
tive states are not equivalent to each other. This reveals 
that L a b 3 is not equivalent to a subfamily of L a bc- 2 ob- 
tained by setting a = c. 



TABLE III. SLOCC classifications of L a b 3 and L a bc 2 









^4 





a = b = 


ab = & a / 6 


a6 / 


Labc-2 


T 2 






(a — c) a = b — 


a=0&6/0 


a/0&6=0 


ab / 



To determine if a four-qubit state belongs to a family 
according to the criteria given by Verstraete et al. [3J 
and Lamata et al. [5], one needs to check if the state 
is equivalent to the representative state of that family. 
For the classification scheme proposed in this Letter, to 
determine if an n-qubit state belongs to a family, one 
needs only to calculate the rank of the coefficient matrix 
of the state. 

In summary, we have studied SLOCC classification for 
general n-qubit states via the invariance of the rank of 
the coefficient matrix and given several examples for n up 
to six. We have also characterized full separable states 
and genuinely entangled states in terms of the rank. We 
expect that the proposed entanglement classification for 
general n-qubit states may find further experimental con- 
sequences. 

This work was supported by NSFC (Grant No. 
10875061) and Tsinghua National Laboratory for Infor- 
mation Science and Technology. 



V. DISCUSSION AND SUMMARY 

Chterental et al. (see Remark 3.5 in [4J) stated that the 
family L a b 3 is equivalent to a subfamily of L a bc 2 obtained 
by setting a — c, where L a b 3 and L a & C2 are given by Q 

L ab3 = a(|0000) + |1111)) + £i_(|0101) + |1010)) 

+ iZ*(|0H0> + |1001)) + 4|(|0001) + 10010) 

+ |0111) + |1011)), (10) 

L ab c 2 = ^-(|0000) + |1111)) + -y-(|0011) + |1100)) 
+c(|0101) + |1010)) + |0110). (11) 

In terms of the rank, the families L a (, 3 and L a i, C2 with 
a = c are both divided into four subfamilies, see Table [TTTT 
As can be seen, the subfamily T<i of L Q {, 3 is a single class 
with representative ^(|Q0Q1) + |0010) + |0111) + |1011», 
whereas the subfamily T% of L a bc 2 ( a = c) is a single class 
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